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We report new results of Bose-Einstein condensation of polaritons in specially de- 
signed microcavities with very high quality factor, on the order of 10 6 , giving the 
polariton lifetimes of the order of 100 ps, which is much longer than their thermal- 
ization time. We see qualitatively different physics from previous experiments with 
polaritons with much shorter lifetime, in part because the longer lifetime in these new 
structures allows the polaritons to travel macroscopic distances up to a millimeter. 
We observe novel effects in three regimes of polariton density. At low density, we 
see classical ballistic motion over hundreds of microns. From this, we obtain a direct 
measurement of the mean free scattering time of the polaritons. At moderate density, 
we see coherent wave-like motion over long distance. This appears as a continuous 
transition, and can be described in terms of quasicondensation. At high density, we 
see a very sharp transition to a trapped state with very high coherence indicative of 
equilibrium BEC. We compare our results to previous experiments on transport of 
polaritons. 



I. INTRODUCTION 

Previous work by many experimental groups has shown that polaritons in microcavities at high 
density and low temperature act in many ways the same as Bose-Einstein condensation (BEC) of 
cold atoms, with numerous resulting effects such as quantized vortices PQ, Josephson oscillations 
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between two condensates with time-dependent period that depends on the particle density [2], 
bimodal momentum distribution OH], indications of a Bogoliubov linear excitation spectrum [SHZ] 
and many other fascinating effects (for reviews see Refs. [5HlU]). Essentially, what is done in these 
experiments is that photons, which are normally massless and noninteracting, are given an effective 
mass by means of a cavity and are given mutual interactions by means of a resonance with an 
electronic excitation in a medium. The effective m these systems reaches world-record values, 
around four orders of magnitude higher than typical nonlinear optical media [11] , and this makes 
the system behave in ways quite different from a standard laser. The mathematical description of 
the interacting photons, that is, the polaritons, is given by the theory of a weakly interacting Bose 
gas in two dimensions, with the addition of generation, drag, and decay terms. 

In the experiments up to now, the lifetime of the particles at low momentum has been of the order 
of a few picoseconds [12|. IT3] , comparable to or just slightly longer than the collision time between 
the particles. This short lifetime does not prevent most of the canonical effects of equilibrium BEC 
from occurring in the polariton gas, because the lifetime can still be longer than the collision time 
between the particles; a number of theoretical |14l IT5] and experimental [3 [16] works have shown 
that a boson system with generation and decay will still have most of the macroscopic effects 
of superfluidity and condensation, and numerical models using a quantum Boltzmann equation 
approach [IT] have shown that a high occupation of the polariton ground state can occur by known 
scattering processes with incoherent, high- kinetic-energy generation of the particles and radiative 
decay in steady state. 

In this paper, we report experiments with new polariton microcavity structures with extraordi- 
narily long-lived polaritons. With lifetimes of the order of 100 ps, more than an order of magnitude 
longer than previously reported lifetimes, the particles can live much longer than their average colli- 
sion time and longer than the time needed to thermalize to the solid lattice temperature by phonon 
emission. Although 100 ps may still seem like a very short time, the absolute time scale is not 
important; what matters in any system is the ratio of the lifetime of particles to the collision time 
between the particles. When the lifetime is much longer than the collision time, equilibrium statis- 
tics apply to a very high degree of approximation. This is the case for atomic condensates |18] , in 
which the atoms have a lifetime in the trap of a few seconds, compared to typical interaction times 
of milliseconds. In principle, the rate of collisions in these polariton experiments can be increased 
by raising the density of the particles, but there is an upper bound on the density of the polaritons: 
their density is limited by a Mott transition at high density to uncorrelated plasma \19\ [20] ; even 
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at lower carrier density than the Mott transition, the reduction of the oscillator strength of the 
exciton due to many-particle screening may decouple the photons from the exciton resonance. In 
the experiments reported here, the lifetime of the particles has been increased while the density and 
collision rate remain approximately the same. These experiments move the polariton system into 
the same range of ratio of lifetime to collision time as the cold atom condensates. The interaction 
between the polaritons is also tunable: it can be made much weaker so that interactions are almost 
turned off; however, it is not yet possible to create an attractive interaction between the polaritons 
or a divergence analogous to a Feshbach resonance. 



II. EXPERIMENTAL METHODS 



The polaritons in these structures exist in a two-dimensional plane, defined by the two parallel 
mirrors which make up the optical cavity and by the parallel sets of GaAs quantum wells between 
the mirrors, placed at the antinodes of the confined optical cavity mode. The only important 
difference between our new structures and the old structures [H EI] is that the cavity formed by 
the distributed Bragg reflectors (DBRs) used to make the microcavities now has a quality factor (Q) 
of around 10 6 , with a designed cavity lifetime of 400 ps. Our previous experiments used a structure 
which was identical in design to this one, but with half as many layers in the DBR mirrors. The 
designed cavity lifetime for that structure was 1.5 ps, with a Q of 4800. The high Q calculated 
for our new structure is confirmed experimentally by the intrinsic line width of the lower polariton 
reflectivity line, which is narrower than our instrumental resolution of 0.05 nm. 

The cavity width in our structures varies continuously as a function of the distance from the 
center of the wafer. Consequently, the cavity photon energy on one end of the structure is below 
the exciton energy in the GaAs quantum wells, while on the other end of the structure, the photon 
energy is above the exciton energy. When the photon energy is below the exciton energy, the lower 
polariton state is mostly photon-like; this is called negative detuning, with the detuning parameter 
5 = -Ephoton — -^exciton- in the reverse situation, with 5 > 0, the polaritons are mostly exciton-like. 
Exactly at the resonance point on the sample where the photon and exciton energies are equal, the 
polaritons are 50-50 quantum superpositions of a photon and an exciton, with an energy splitting 
between two polariton states known as the Rabi splitting, which in this case is approximately 13 
meV. In our samples this Rabi splitting is large compared to the intrinsic line widths of the photon 
and exciton modes; this is known as the strong coupling limit. 
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The interaction of the polaritons with each other and with the phonons in the lattice comes 
entirely through their excitonic part. Therefore, by moving to a region of the sample with negative 
detuning, we can reduce the interactions. There is a tradeoff in what detuning to choose. When 
the particles are more excitonic, they will thermalize much better to the lattice; indeed, they will 
reach the lattice temperature. Polaritons reaching nearly the lattice temperature when they are 
excitonic has already been reported for short-lifetime polaritons [22]. In this regime, however, the 
mean free path, and therefore the diffusion length, of the polaritons will also be shorter. For the 
experiments discussed here, we have chosen a region on the photonic side of the sample where the 
polaritons have relatively weak but nonzero interactions with the lattice and with each other. This 
allows them to have long distance motion. 

In all of the experiments discussed here, the excitation laser photon energy is well above the 
polariton energy, with wavelength 707 nm, at the third minimum in reflectivity above the micro- 
cavity stop band (see Appendix A for the reflectivity spectrum). Since the excess energy is large, 
many phonons were emitted by the excited carriers before they turned into polaritons; as a result, 
the coherence of the pump laser is lost to the incoherent phonon bath and not transmitted to 
the polaritons. The bath temperature in which the sample was held was approximately 7 K. A 
Princeton intensified CCD camera on a 0.3 m spectrometer was used for time-integrated spectral 
images; a Hamamatsu streak camera was used for time-resolved luminescence. 



III. LOW-DENSITY BEHAVIOR 

Figure [T] shows data for low-density polaritons created at a spot on the sample on the photonic 
side of the resonance (detuning 5 = —10 meV), with very low excitation density. The shift of 
the lower polariton energy with detuning due to the gradient in the cavity width gives a spatial 
gradient of the potential energy felt by the polaritons. The spatial motion of the polaritons seen in 
Fig. [ijb) is at first surprising: the polaritons move uphill, to higher potential energy, for hundreds 
of microns; uphill motion more than 1 mm has been observed in this structure. This uphill motion 
can be understood if we take into account that the laser which generates the polaritons is tuned to 
a photon energy well above the highest-energy polariton states. As discussed above, free electrons 
and holes are generated at very high kinetic energy and then lose energy by phonon emission, 
eventually entering the exciton and polariton states. Therefore polariton states with a broad range 
of momenta are occupied near the laser excitation spot. Some of these polaritons drop down to lower 
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energy by phonon emission and polariton-polariton scattering, but many propagate outward from 
the laser excitation spot without much scattering. Those with momenta pointing in the direction 
of the uphill gradient of potential energy will trade their kinetic energy for higher potential energy 
until they eventually hit a place where the potential energy of the polaritons equals their initial 
energy, as illustrated in Fig. [ljd). At this point of zero momentum, they can turn into photons 
which are emitted normal to the cavity plane. These photons are what we record for the data 
shown in Fig. [ljb) . 

If we resolve the photoluminescence by angle of emission, we get a fuller picture of what is 
happening. Fig. [ljc) shows the emission spectrum resolved by angle of emission for very similar 
conditions to Fig. [ljb) . The angle of emission has a one-to-one mapping to the in-plane momentum 
of the polaritons before they turn into external photons. As seen in this figure, the emission 
corresponds to a parabola which is smeared to the left. The initial parabola (indicated by the 
white line) is the dispersion E{k) of the polaritons at the location of the laser spot, calculated 
by a fit to the complete set of lower polariton dispersion data across the sample. This dispersion 
gives the effective mass of the particles according to E = Eq + (hk\\) 2 /2m, where h\ is the in-plane 
momentum of the polaritons and Eq is the low-density, zero-momentum energy of the polariton. 

The smearing comes from the force on the polaritons due to the cavity gradient, which causes a 
shift in the polariton energy. Since F = —dU/dx and F = Hdk/dt, the particles feel an acceleration 
to lower potential energy and to negative k. Particles starting out in the uphill direction slow down 
but predominantly remain at the same energy until they reach k = 0, because they have weak 
scattering and therefore move ballistically. Some even continue past k = 0, turning around to go 
back the other way. Particles which start in the downhill direction, with negative k, accelerate in 
the same direction and move to ever larger magnitude of k, and so are never seen in our real-space 
imaging with k ~ 0. 

Figure [2^a) shows time-resolved data of the same type of motion. As seen in this figure, there 
is a time delay for the light emission coming from locations on the sample far away from the laser 
spot. This occurs because the polaritons take time to propagate from the creation spot to the 
observation spot. Fig. [2](b) shows the time delay for the onset of the light emission as a function 
of the spatial separation from the generation spot. The theory curve is the predicted transit time 
from a semiclassical model, using the known force (equal to 8.5 meV/mm at the excitation spot and 
decreasing in the direction towards the zero-detuning location on the sample) and the dispersion 
E{k) of the polaritons. The fact that the polaritons do not have constant mass at large k must be 
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taken into account to generate the fit for long propagation distances (see Appendix B). 

The data of Fig. [2^a) are well fit by a single exponential decay convolved with a Gaussian time 
spread of 75 ± 10 ps. This time spread increases only slightly for the time-resolved data far from the 
laser generation spot. This indicates that the spread in time is almost entirely due to the spectral 
diffusion which occurs at the laser spot as the carriers lose the initial excess energy given to them 
by the laser photons (about 120 meV) as they fall toward the polariton states. Very little additional 
time spread occurs as the polaritons move hundreds of microns across the sample, consistent with 
the view that they are moving nearly ballistically. 

Fig. [2^c) shows the lifetime of polaritons traveling ballistically at constant energy, using the 
data of Fig. Qc) for the intensity as a function of momentum and the time-of-fhght data to map a 
momentum change to a time-of- flight. At each energy, the data fit a single exponential. As seen in 
Fig. [2^c), the effective decay time decreases with increasing polariton energy. This effective decay 
time measures the total rate of departure from the ballistic states, i.e., 1/t c s = 1/tr + l/r sca tt, 
where tr is the radiative lifetime and r sca tt is the average scattering time for all scattering processes 
(phonons, defects, and other polaritons) which enters into the mean free path, I = ur sca tt- A full 
quantitative analysis of the energy dependence of the scattering is beyond the scope of this paper, 
but part of the variation with energy (about 20%) can be understood as due to the increase 
in the excitonic fraction of the polaritons with increasing energy, which causes them to interact 
more strongly. Both polariton-polariton scattering and elastic polariton scattering with defects are 
independent of the energy of the polaritons in 2D; polariton-acoustic phonon scattering, which is 
proportional to E 2 for emission of phonons into the 3D bulk [T7], should be very weak. 



IV. INTERMEDIATE DENSITY BEHAVIOR 

Figure [3] shows data for polariton motion under similar conditions as Fig. [TJ but at higher laser 
pump intensity, which corresponds to higher polariton particle density. As seen in the momentum- 
space data shown in Fig. [3^c), the energy distribution of the particles has become compressed to a 
single energy. This energy is blue-shifted relative to the original low-density ground-state value. 

The blue shift occurs because the off-resonant laser pump creates not only polaritons but exci- 
tons, and this cloud of excitons has an overall repulsive interaction with the polaritons which leads 
to an additional potential energy felt by the polaritons. Since the excitons have mass of the order 
of an electron mass, about 10 times larger than the mass of the polaritons, the excitons are essen- 
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tially static as seen by the polaritons, and form a fixed potential energy barrier for the polaritons. 
The dashed line in Fig. [3^b) shows the spatial potential energy profile felt by the polaritons in the 
presence of the exciton cloud and the gradient in the cavity width. The polaritons generated by the 
laser begin their existence sitting on top of the nearly-static exciton cloud. In the data of Fig. [TJ 
the potential energy due to the exciton cloud was negligible, but at the higher excitation density 
of Fig. [3j the potential energy shift of 1.9 meV due to the exciton cloud is much larger. 

After they come into existence sitting on top of the exciton cloud, the polaritons can accelerate 
away from this spot. This motion does not come from the broad range of high energies and 
momenta seen in Fig. [T] at low density, because the polaritons in this higher-density regime now 
start their motion almost exclusively in the ground state at the excitation spot, with nearly zero 
initial momentum at that location. The highly peaked potential then causes them to accelerate 
both uphill and downhill. The uphill and downhill motion correspond to positive and negative 
peaks in momentum space, as seen in Fig. (3^c). The particles initially accelerated in the uphill 
direction will slow down until they hit the point in the sample which has the same potential energy 
as their starting point, at which point they will have zero momentum, and then they will turn 
around to accelerate back downhill. This point at which they turn around can be clearly seen in 
the spatially-resolved data of Fig. [3^b). The total energy of the particles does not change during 
this exchange of potential and kinetic energy, as expected when there is no drag. 

Although the data of Fig. [T] and Fig.[3]look quite different, the spatial and momentum profiles at 
low density evolve continuously from those of the type shown in Fig. [T] into the type seen in Fig. [3] 
as the particle density is increased. As the density is increased, the polaritons are increasingly 
attracted to a single, spatially extended, energy state resonant with the potential energy at the 
point of creation. Eventually, this energy state dominates the entire spectrum. 

We can understand this as arising from the boson nature of the polaritons, which causes them 
to be attracted into the same state as the occupation number increases. It is worth some discussion 
whether one can consider the state of Fig. [3] to be a superfluid state. We can say, first of all, that 
the motion occurs without drag or damping, since the energy of the particles does not change; it is 
constant over the full range of motion over 200 /xm, as seen in Fig. ^h). In addition, we can say that 
this motion is coherent. This is seen first of all in the spectral narrowing, which indicates increased 
coherence. This width is controlled by small fluctuations of the laser intensity, which lead to changes 
in the exciton cloud potential. Interference measurements using a Michelson interferometer also 
show that the light emitted from the creation spot is phase coherent with the light emitted from 
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the turnaround point, more than 200 //m distant. Fig. j4|a) shows that the coherence measurement 
gives a coherence time of about 40 ps, but the maximum coherence is offset by 140 ps, which 
corresponds to the propagation time from the creation spot to the turnaround point. This is a 
fascinating example of propagated coherence. 

A momentum-space distribution similar to that of Fig. [3^c) has been reported by the Grenoble 
group [23\ I24|. under very similar conditions of non-resonant excitation, but with short-lifetime po- 
laritons. The effect occurred for the same reason, namely, the acceleration of a coherent superfluid- 
like state away from the potential energy peak created by the exciton cloud. In our experiments, the 
momentum distribution is dominated by the fast-moving particles leaving the excitation spot; so 
the long lifetime of the polaritons does not make a strong difference in the momentum distribution. 

The long-distance flow in real space in our experiments is strongly affected by the lifetime, 
however. Fig. [4](b) shows fits of the spatial distribution to a solution of the one-dimensional 
Gross-Pitaevskii equation for a coherent superfluid with negligible particle-particle interaction, 
with generation at the laser spot and an overall lifetime decay, for the potential energy profile given 
by the exciton cloud and cavity gradient. The details of this model are discussed in Appendix B. As 
seen in this figure, the fit is quite good when a lifetime of 113 ps is used. The fit is highly sensitive 
to the lifetime of the polaritons, since they must live long enough to travel across the sample to 
the turnaround point. This lifetime is also an underestimate of the true lifetime, because flow in 
the sideways direction, which occurs because the system is actually two-dimensional, will give an 
effectively shorter lifetime for the flow in the uphill direction. 

Most of the canonical effects of superfluidity, such as quantized vortices and Josephson oscil- 
lations, are the result entirely of long-range phase coherence, which allows the superfluid to be 
modeled by a Gross-Pitaevskii equation, as seen here. Many would add an additional criterion 
for the proof of superfluidity, namely a critical velocity at which the superfluidity (i.e., long-range 
phase coherence) ceases. We have not measured a critical velocity. It is clear, however, that the 
phase coherent state is stable and resistant to dephasing and drag, as evidenced by the fact that 
the phase coherence is maintained as the fluid flows hundreds of microns at constant energy. 

Long-range coherent transport of polaritons has been seen in other experiments. In the ex- 
periments of the Madrid group [25], short-lifetime polaritons were seen to move with soliton-like 
behavior over long distances. Although the lifetime of the polaritons in those experiments was very 
short, the polariton population was continuously replenished by a large laser spot which blanketed 
the entire space where the polaritons traveled. In our experiments reported here, the creation 
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of polaritons occurs only at the laser excitation spot, and there is no additional pumping as the 
polaritons travel. 

The group in Paris [22] has shown coherent motion of polaritons of around 50 fim in quantum 
wires, with polaritons with lifetime of around 15 ps. In those samples, the scattering of the po- 
laritons was suppressed due to the one-dimensional motion. In our experiments reported here, the 
polaritons move in a two-dimensional plane, with very low scattering due to the high quality of 
our structures. As in our experiments reported here, the Paris group used the nearly static exciton 
cloud as a potential barrier for the polaritons; this method of creating a potential barrier has also 
been demonstrated by the Dortmund group [27J and the Cambridge group [28J with short-lifetime 
polaritons, and is a very promising method for control of the motion of polaritons. 

Incoherent motion of polaritons over tens of microns has also been reported for short lifetime 
structures. In Ref. [29J, downhill drift of up to 50 /im was seen in a cavity thickness gradient. 
In Ref. [1], drift of up to 40 /xm was seen for a force of approximately 50 meV/mm created by 
a strain-induced shift of the exciton energy. This was enough to cause the polariton distribution 
created by the laser to shift to have a local maximum at the center of a harmonic potential trap, 
away from the center of the laser peak. In those experiments, in the absence of the drift force, the 
diffusion distance was less than 10 /im. 

V. HIGH DENSITY BEHAVIOR 

Figure [5] shows data for polariton spatial and momentum distributions under the same conditions 
as Fig. [3j but at even higher density. Again, there is a very dramatic change. The energy of the 
polaritons decreases, and the momentum distribution and spatial distribution are very compact. 
The product of the widths of the momentum and space distributions is 1.8 times the uncertainty 
limit. The total intensity of the light from this spot jumps nonlinear ly, and the emitted light 
becomes coherent. The spectral width of the light in this case is limited by the resolution of 
our detection system, at 0.05 nm. The coherence time of this spot, measured by a Michelson 
interferometer, is at least 280 ps; this number is limited by fluctuations of the laser power, although 
it is less sensitive than the intermediate-density case of Fig. [3j because it does not sit on top of the 
exciton cloud. 

In this case, the polaritons are now spatially trapped. Although the trap profile which we create 
is not a confining potential such as a ring [30j , the combination of the exciton cloud and the sloped 
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cavity gradient still appears to confine the polaritons in both dimensions at the highest density, 
since in spatial imaging measurements, the spatial profile of the state seen in Fig. [5^b) is localized 
in both spatial directions to a radius of approximately 5.5 /mi. As seen in this figure, the bright 
emission at high density does not come from the laser excitation spot; it comes from a region about 
16 /im to the right, at the energy minimum of the trap created by the wedge of the cavity and the 
exciton cloud. 

The jump to a trapped state may be related to one of two mechanisms. One possibility is that 
the initial velocity given to the coherent polaritons exceeds a critical velocity, and they can no longer 
flow long distances, and therefore fall into the potential energy minimum where they are trapped. 
Another possibility is that the threshold for the trapped state is the density at which the particles 
undergo "true" Bose-Einstein condensation, and the long-range coherent behavior seen at lower 
density is a "quasicondensate" with increasing density. Although there is a discrete jump down in 
the energy of the polariton emission from Fig. [3^b) to Fig. [5^b), this jump does not correspond to 
a jump between discrete energy states in the trap. The size of the trapping region shown in Fig. [5] 
gives states separated by around 30 /ueV, i.e., nearly a continuum of states when compared to the 
typical energy scale of meV in these data. The jump is between a state resonant with the creation 
energy of the polaritons on top of the exciton cloud, as seen in the data of Fig. [3j and the ground 
state of the trap. 

The fact that the energy of this emission is lower than the energy of the emission at lower pump 
powers shows that this emission is not standard lasing — standard lasing has been shown to occur in 
this type of structure, and always occurs at higher energy, at the cavity photon mode [2T | l3T | 132]. 
In the data shown in Figs. [3] and [5j the bare cavity photon energy is at 1.60 eV, well above the 
relevant polariton energies. Also, the fact that the trapped state is spatially displaced from the 
location of the laser spot by about 16 /im, and is in a region with almost no pump laser intensity, 
also indicates that standard pumped lasing is not occurring. 

We interpret this highest-density regime as "true" Bose condensation in a trapping potential. 
Although true Bose condensation cannot occur in a translationally invariant two-dimensional sys- 
tem, it can occur in two dimensions when there is a trapping potential [63\ I34| . The combination of 
the exciton cloud and the gradient of the cavity energy makes a trap. At a critical density threshold, 
the long-range superfluid-like motion seen in Fig. [3] is superseded by flow into the trapping region. 
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VI. THREE REGIMES 

Figure [6] shows the density dependence of several parameters for conditions like those of Figs.[TJ 
[3j and [5] In Fig. [(^a), we show the total intensity of the ground state luminescence as a function 
of excitation density; this plot gives the total intensity at k\\ ~ 0, taken from momentum-space 
data like that shown in Figs. [T](c), [3^c), and[5](c). In Fig. [6^b), we show the fraction of the total 
luminescence which is emitted from a narrow spectral range near the ground state, integrated 
over all momenta. In Fig. [6^c) we show the spectral width (full width at half maximum) of the 
luminescence. 

The three regimes discussed above have distinctly different behavior. In the lowest density 
range, the luminescence intensity increases linearly with the excitation intensity. The fraction of 
the particles at low energy remains constant, equal to approximately 0.1. In this regime, the energy 
distribution of the polaritons is controlled by their emission of phonons, which is weak since in this 
region of the sample they are mostly photon-like. 

In the middle regime, the fraction of particles in the ground state increases superlinearly with 
excitation density. The ground state is quasi-coherent, and increasingly has long-range flow in a 
single energy state. We interpret this regime as occurring when collisions between the particles 
become significant, so that the polariton gas becomes more thermalized, and boson stimulated 
scattering favors scattering into low-energy states, producing a quasicondensate which increases 
continuously with density. 

The upper density regime occurs very abruptly, and corresponds to an order of magnitude 
greater intensity in a much more compact spatial and spectral distribution. This abrupt behavior 
is characteristic of a true phase transition. We believe that the transition is in fact more abrupt 
than seen in these data; although we use a diode laser with intensity stabilization, the excitation 
intensity still has some fluctuations, of the order of the width of the transition seen here. When 
a stabilized laser is used, the spectral width, and consequently the measured coherence time, is 
greatly increased [35j . 

Several previous experiments O l3T| [361 [37] with short-lifetime polaritons have shown a more 
gradual transition. In experiments by the Paris group [3B]> the sharpness of the transition could 
be varied by changing the volume of the sample, similar to how the S-curve of a laser onset can 
be sharpened by changing the volume [391 140] . We emphasize that in the present experiments, 
the sharpness of the transition we see at high density is not related to a volume change, since the 
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volume is exactly the same as that of previous structures [4] which had a much broader transition 
range; only the quality factor of the cavity is changed. We believe that the increased sharpness of 
the transition seen here is connected to the better equilibration of the particles due to their long 
lifetime. 

VII. FUTURE DIRECTIONS 

Because the trap in these experiments comes about due to the combination of the cavity gradient 
and the exciton cloud, it is difficult to model the dynamics of the transition from ballistic classical 
motion to coherent long-range motion to a trapped state with any exact theory. Therefore we are 
proceeding to do experiments with a controlled harmonic potential trap created by inhomogenous 
stress, as in previous experiments [U 01], and we will report the results of the experiments using 
a controlled harmonic potential in a future publication. Figure [7] shows that when we use a stress 
trap, we can make two trapped condensates side by side. In this case the exciton cloud acts as a 
barrier between the two sides of the stress trap. 

We have also seen evidence of spontaneous vortices in this system, which will also be reported 
elsewhere. In general, there are many exciting new directions to take with the combination of 
long lifetime polaritons, the stress trap, and the exciton cloud as a movable potential barrier. It 
is also possible now to use the AC Stark effect to create a potential for the excitons on very short 
time scales |42| . The new structures with long lifetime open up qualitatively new physics on a 
macroscopic scale. 

Acknowledgements. This material is based upon work supported by the National Science 
Foundation under grant DMR- 11 04383. The work at Princeton was partially funded by the Gordon 
and Betty Moore Foundation as well as the National Science Foundation MRSEC Program through 
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Appendix A: Additional Experimental Details 

Fabrication. The high-Q structures used for these experiments was identical in design to that 
used in Ref. [4| except that the number of layers in each of the distributed Bragg reflectors (DBRs) 
was doubled to 40 layers on the substrate side and 32 layers on the top side. This required over 
30 hours of molecular beam epitaxy (MBE) time. The rate of MBE growth during this time was 
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established to be constant to within 2% by comparison of the growth parameters before and after 
the fabrication and by comparison to optical reflectivity data. 

Excitation. For incoherent excitation, we use excitation wavelength at the third dip above the 
reflectivity stop band of the cavity. This generates hot electrons and holes which then fall into the 
exciton and polariton states by phonon emission (see Ref. [T71 ) 

We used a stabilized diode laser to avoid fluctuations of the pump intensity; these fluctuations 
have been shown |35j to lead to significant fluctuations in the spectral line position of the polaritons 
and broadening of their spectral line width. The tradeoff of using a stabilized diode laser is that it 
had maximum power of 50 mW, which limited the highest power to that shown in Fig. 6. 

Momentum-space measurement. The far- field angular distribution of the photolumines- 
cence from the sample gives the momentum distribution of the polaritons, since the generalized 
SnelFs law requires that the in-plane momentum of the polaritons be conserved when they convert 
to photons outside the system, and this in-plane momentum maps to the angle of photon emission 
in the Fraunhofer far- field limit. We used a lens system to project the far- field emission onto the 
entrance slit of our Acton 0.3-meter spectrometer, which has a Princeton Instruments intensified 
CCD camera. This allows either energy vs. momentum images as direct data in first order transmis- 
sion through the spectrometer, or 2D momentum-space images in zero-order transmission through 
the spectrometer. This system is very fast since all momentum-space data is taken in parallel by 
the CCD camera. 

Calibration of the fc-space emission data. The k = point of the angle-resolved lumines- 
cence measurements was found by bouncing a laser from the surface of the sample to find the angle 
which gave exact retroreflection. This method had an uncertainty of ±10 3 cm -1 . The E(k) dis- 
persion of the polaritons at low density is known from the Rabi splitting; this was checked in many 
previous calibrations of /c-space data. The mapping of the measured angle to k\\ in cm -1 could 
therefore be done using a fit of this dispersion relation to the low-density angle-resolved data. This 
mapping was also checked by a physical measurement of the angle of the rays using an adjustable 
iris. 

Spatial selection of time-of- flight data. For the data of Fig. 2(a), an imaging system 
was used to project the real-space distribution of the photon emission onto the entrance slit of 
the spectrometer with a streak camera. The different curves of Fig. 2(a) correspond to different 
wavelength selections of the spectrometer. Since, as seen in Fig. 1(b), the energy of the k\\ — 
emission maps to the position in the uphill motion of the polaritons, this method selects out the 
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emission at different points on the sample. 



Appendix B: Additional Details of Analysis 

Theoretical model of polariton motion. The one-dimensional Gross-Pitaevskii equation 
with decay is 

.. di\) ft 2 d 2 ip T ,. , , ,,, ,,2 , ih , 

where ^ is the macroscopic wave function of the superfiuid, V(x) is the external static potential, U 
is the particle-particle interaction strength, and r is the particle lifetime. In the low-density limit, 
the nonlinear term becomes negligible, and this equation becomes simple the Schrodinger equation 
for a macroscopic wave function. The external potential included the gradient due to the the cavity 
plus the exciton cloud potential, which for the data of Fig. 3 was a Gaussian with a full width at 
half maximum of 14 ^m and height of 2 meV. To solve this equation, we determined the eigenstates 
of a polariton on the potential profile taken from the experimental parameters, consisting of a slope 
plus a Gaussian peak by numerically diagonalizing the spatial Hamiltonian. An initial Gaussian 
wavepacket was created at the peak of the perturbation, with an energy width of roughly 0.15 
meV. The wavepacket was evolved according to its constituent eigenstate. To simulate steady 
state conditions, the real space and k-space distribution of the wavepacket was integrated over the 
evolution time of the simulation, which was 2.8 times the lifetime. 

Average lifetime. The average lifetime of the polaritons in a Maxwell-Boltzmann distribution 
is given by 



- = — [ e-( E - E °y kBT C 2 (E)D(E)dE, 

T T Dhot J 



Tphot 

where r p hot is the cavity photon lifetime, C(E) is the Hopfield coefficient giving the fraction of 
the polariton state which is photonic, and D(E) is the density of polariton states at energy E. 
C(E) and D(E) are calculated knowing the Rabi splitting of the polaritons and the cavity photon 
dispersion relation, assuming the excitons have constant energy. This calculation assumes that the 
excitons have negligible decay rate, which is valid since the excitons decay primarily by scattering 
down into polariton states, with negligible nonradiative decay in these samples. 

In Ref. HI the cavity photon lifetime was 2 ps, which gave an average polariton lifetime of 16 
ps at the effective temperature of the polaritons of 30 K. The average lifetime is longer than the 
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cavity lifetime because the radiative rate of higher-/c states drops due to their decreased photonic 
component, when the k = state is at resonance. 

Semiclassical model for ballistic motion. The low-density data of Fig. 2 in the text was 
fit to a model in which we assumed classical motion of particles but with E(k) varying according 
to the calculated lower polariton spectrum instead of a single mass. The formula for E(k) for the 
lower polariton is 



ELp{k) — -(-Ecxciton + -Ephoton) _ 2 V ^ cxciton ~~ -^photon) 2 + 4^ 2 , 

where Q is the Rabi frequency, equal to 1/2 the gap between the upper and lower polariton modes 
at resonance. The exciton energy -Ecxciton is assume constant over the relevant momentum range, 



effective cavity length, which depends on the position on the sample, namely the distance r from 
the center of the wafer. 

For the fixed-mass model, one can solve for the time of flight according to Newton's equations. 
F = — W(x) = ma = m(d 2 x / 'dt 2 ) , where V{x) indicates the k = energy of the polariton across 
the sample due to changing cavity photon energy but can be though of as a changing potential 
energy. This is simple to apply if we approximate the slope as constant. The polaritons are starting 
with initial momentum and rolling uphill until they lose all kinetic energy and are observed in their 
k=0 state. In this way we can easily find the time of flight according to the mass, final energy 
and initial k=0 energy of the polariton. This method fails because the mass of the polariton is 
not constant during its motion. The effective mass picture breaks down at large momenta where 
the polariton dispersion relation turns over from being steep and photon-like to being flat and 
exciton-like. For these reasons, we expect this model to break down at larger distances. 

If we focus on momentum, F = —W(x) = dp/dt, based on knowledge of the initial momentum 
we can calculate the time of flight until the momentum equals zero. This means that we can 
directly assign a time of flight to each energy and therefore distance. This method is only trivial 
when VV(x) = constant. Otherwise we are forced to numerically integrate iteratively between 
steps in x and p, since the velocity depends on p. 

The next logical step is to apply the full Hamilton's method of classical dynamics. We can 
directly write down the full Hamiltonian of the polariton based on knowledge of its parameters 
from other experiments. The relevant model is as follows: 



1 




and the photon energy is given by -Ephoton 




where L(r) is the 



E x (x,p) = E + ax + 




■x 
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2 

Ep(x,p) = E + bx + 



nEo+bx 

Z -2/r,2 



c 2 /" 2 

„ , s E x (x,p) + Ep(x,p) ^(E x (x,p)-Ep(x,p)y + 4W 

ttLP(X,p) = • 

Once we have the polariton parameters (deduced by observing the PL energy versus position on the 
sample) we can use Hamilton's equations of motion to numerically integrate for the time of flight. 
This method fully accounts for the changing slope of the polariton as well as the changing effective 
mass. This method completely eschews the concept of effective mass and properly accounts for the 
non-parabolic dispersion relation at moderate-high momenta. 

The polariton parameters that yielded the best fit to the time-resolved data were En = 1.6055 eV, 
f2 = 5 meV, a = 2 meV/mm, b = 14 meV/mm. These are slightly different from the parameters 
deduced from a direct fit of the PL energy versus energy, which gave = 6.5 meV, within the 
uncertainty of the parameters for this sample. 

Time t = was defined as the limit of the fit value of to for the curve e~ <yt ~ t °^ T as x — > 0, for 
fits to the data like those shown in Fig. 2(a). 
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FIG. 1: Long-distance ballistic motion of polaritons at low density, a) Spatially-resolved photon emission 
from hot carriers created by the laser excitation spot hitting the microcavity structure, at a spot on the 
sample with 5 = — 10 meV (the polaritons are mostly photon- like) , in the low-density limit (1 mW pump 
laser in a spot of radius 12 fim). This shows the spatial extent of the excitation region where the polaritons 
are generated, b) Spatially-resolved photon spectrum from the lower polaritons under the same excitation 
conditions. Only photons emitted with angle ±2.5° relative to the normal to the surface are detected, 
corresponding to in-plane polariton momentum ku ~ 0. The polaritons generated at the excitation spot 
move away in all directions, but in the fcy ~ emission we see only those which move ballistically in the 
uphill direction, slowing from high momentum to low momentum, c) Momentum-resolved spectrum of the 
photon emission from the lower polaritons under the same excitation conditions but slightly less detuning 
from resonance (5 = —3.2 meV; the pump excitation occurs at the spot where the lower polaritons have 
energy 1.596 eV). The smearing of the distribution toward the right corresponds to the momentum gained 
by the particles as they accelerate downhill. The real-space distribution for detuning S = —3.2 meV is given 
in the supporting information, d) Illustration of the process by which the polaritons appear to move uphill. 
Polaritons in high-fc states which escape the excitation region travel ballistically. As the cavity width shrinks, 
the potential energy of the polaritons increases, and their kinetic energy drops, until they reach k\\ = 0. 

FIG. 2: Time- of -flight data for the long-range ballistic motion, a) Blue lines: Time-resolved intensity of 
polariton emission for three spots on the sample at different distances from the excitation spot, following a 
short (~ 5 ps) laser pulse at 725 nm. The excitation conditions were similar to those of Fig. [ija) (20 mW 
peak power in a spot of radius 14 /nm, at a location with lower polariton energy equal to 1.5857 eV.) Black 
lines: hot luminescence signal (which follows the laser pulse) giving the instrumental time resolution. Red 
lines: fits to a simple theory of a a single exponential decay e~^ t ~ t °^ r convolved with a Gaussian e - * 2 /"' 2 
broadening. The offset times to correspond to the time-of-flight for the polaritons from the excitation spot 
to the spot where the photon emission is detected, b) Black symbols: The time-of-flight, from data like 
that in (a), as a function of the distance from the excitation spot. The three theory curves correspond to 
different levels of approximation for the effective mass and local potential gradient felt by the polaritons. 
c) The lifetime of the polaritons in ballistic states as a function of energy, mapping I(k) from equal-energy 
slices of the data of Fig. [ljc) to I(t) using k(t) from the Hamilton's approach theory developed to fit the 
data of (b). Inset: typical fit of a lifetime to the data, for polaritons at energy 1.5976 eV. 
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FIG. 3: Long-distance coherent motion at high density, a) Photon emission from hot carriers created by the 
laser excitation spot on the microcavity structure, as in Fig. [lja). The excitation conditions are the same 
as for Fig. |ljc) except the pump power is 34 mW. b) Spatially-resolved photon emission spectrum, using 
the same method as Fig. |ljb) , at this higher excitation density. As seen in this figure, there is a blue shift 
of the polariton luminescence at the excitation spot due to the compact, slow-moving exciton cloud at that 
location, which repels the polaritons. At a point further away uphill, which is resonant in energy with the 
potential energy at the top of this exciton cloud, the polaritons which have moved away slow down to ku ~ 
and become visible as a bright spot at the turnaround point. The dashed white line shows the potential 
energy felt by the polaritons, a combination of the gradient due to the cavity width variation and the exciton 
cloud potential, c) Momentum-resolved photon emission from the polaritons under identical conditions as 
the data of (b). The low-energy tail just barely visible in Fig.Qb) has become dominant, giving a horizontal 
line which corresponds to motion at constant energy across the sample. 



FIG. 4: Analysis of the coherent motion, a) Dots: The fringe visibility as a function of time delay for 
interference of the light from the source spot with the light from the secondary spot 200 ^m distant, for 
the two spots shown in Fig. |3jb). Solid line: fit to a Gaussian with half width at half maximum of 37 
ps. The peak coherence is for a delay of 140 ps, corresponding to the propagation time for the polaritons 
between the two spatial spots. There is strong scatter of the data points because of fluctuations in the 
spatial position of the second spot due to pump laser intensity fluctuations. Inset: a typical interference 
pattern, b) Black dots: the intensity as a function of position, from the spatially-resolved data of Fig. |3tb) 




integrated over energy. Colored curves: steady-state spatial distribution from a one-dimensional solution 
of the non-interacting Gross-Pitaevskii equation (macroscopic-wave Schrodinger equation) for the potential 
energy profile felt by the polaritons, given by the dashed white line in Fig. |3](b), with generation only at the 
excitation spot, for various lifetimes. 



FIG. 5: Trapped condensate above a critical density. a)-c) show the hot-carrier luminescence (excitation 
spot), spatially-resolved spectrum and momentum-resolved spectrum under identical conditions as Fig.^a)- 
(c), but with high power (44 mW). As seen in (b), the spatial distribution is now very compact and is located 
at the potential energy minimum of the potential profile (indicated by the dashed white line) . As seen in 
(c), the momentum distribution is also compact. The product of the real space and momentum space full 
widths is 1.8 times the uncertainty limit. 
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FIG. 6: Density dependence in the three density regimes, a) Intensity of the photon emission peak at = 
(emission normal to the surface of the sample) as a function of excitation power, for conditions very similar 
to those of Figs. T][5 (Laser spot size of 8 /im diameter, excitation spot at the point on the sample with 
lower polariton energy = 1.5954 eV.) Thin dashed line: linear dependence (slope = 1 on the log-log plot), b) 
The total fraction of luminescence contained in a spectral range of 0.5 meV centered on the wavelength of 
the peak emission intensity, integrated over all momenta, for the same conditions as those of (a) . c) The full 
width at half maximum of the spectral peak, from the same spectra as used for (a). The spectral resolution 
was 0.1 meV. The vertical dashed line at low density in each case indicates the point of deviation from 
linearity, from the data in (a); the vertical dashed line at high density indicates the threshold pump power 
for the compact condensation shown in Fig. [3J 



FIG. 7: Initial results for two condensates in adjacent traps, a) Spatially-resolved spectrum for conditions 
similar to those of Fig. f5fb), but with an external trapping potential created by a stressor, b) Estimated 
external trapping potential (blue curve) due to the applied stress, and the exciton cloud potential (green 
curve) for this experiment. 
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